Time-dependent density-functional theory with self-interaction correction 
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We discuss an extension of time-dependent density-functional theory by a self-interaction correc- 
tion (SIC). A strictly variational formulation is given taking care of the necessary constraints. A 
manageable and transparent propagation scheme using two sets of wavefunctions is proposed and 
applied to laser excitation with subsequent ionization of a dimer molecule. 
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Density Functional Theory {Y)YT)L-2^ has 
evolved over the last decades to a standard the- 
oretical tool for the description of electronic prop- 
erties in many physical and chemical systems, es- 
pecially in systems with sizable numbers of elec- 
trons. The extension of DFT to Timc-Dcpcndcnt 
situations (TDDFT)^i^ii is a more recent achieve- 
ment still motivating many investigations, both for 
formal and practical aspects^. It turns out that 
TDDFT is one of the few, well founded theories, 
allowing to describe dynamical scenarios in com- 
plex systems, which is a key issue for understand- 
ing microscopic mechanisms, beyond mere ener- 
getic considerations. This is especially true if elec- 
tron emission comes into play as, e.g., in case of 
irradiation processes. 

A practical treatment of DFT, and even more so 
TDDFT, requires simple approximations to the 
exchange and correlation functional. The simplest 
one is the Local Density Approximation (LDA) 
which proved very useful in calculations of struc- 
ture and low-amplitude excitations (optical re- 
sponse, direct one-photon processes)^. It is also 
often used as a first order approach for more vi- 
olent dynamical processes where huge energy de- 
posits lead to a large number of emitted electrons 
as, e.g., in clusters subject to intense laser fields 
or collided by highly charged particles^. How- 
ever, LDA is plagued by a self-interaction error be- 
cause its Kohn-Sham field involves the total den- 
sity including the particle on which the field just 
acts. As a consequence, the LDA produces the 
wrong Coulomb asymptotic of the mean field and 
thus underestimates the ionization potential (IP) 
of a system. This spoils, e.g., the dynamical de- 
scription of excitations involving ionization pro- 
cesses, in particular close to the threshold. A cor- 
rect treatment requires a self-interaction correc- 



tion (SIC). Such a SIC complementing LDA for 
static calculations was proposed ivrt^ and has been 
used since then at various levels of refinement for 
structure calculations in atomic, molecular, clus- 
ter and solid state physics, see e.g .^^d3ii4,i5 _ 

The original SIC scheme leads to an orbital depen- 
dent mean field which causes several formal and 
technical difficulties. There are attempts to cir- 
cumvent these problems by treating SIC with op- 
timized effective potentials (OEP), se o^^d"^ . That, 
however, overrules some crucial physical features 
of SIC, particularly the trend to produce localized 
single-particle states^^. 

Application of SIC in time-dependent situa- 
tions are mostly done in approximate manner, 
linearizedi^, using averaged-density SIC^, or 
relying on various versions of time dependent 
OEP-KLlS2i2i^. rpj^g TDOEP-KLI, however, 
also suffers from inconsistencies as it leads, e.g. 
to violation of zero force theorem and energy 
conservations^. The aim of this letter is to present 
a thorough variational formulation of fully fledged 
TDSIC without further approximation, together 
with a manageable propagation scheme which al- 
lows to obey all boundary conditions, namely the 
zero-force theorem, conservation of energy and or- 
thonormality of the occupied single-particle or- 
bitals. This will in particular serve as a bench- 
mark for the development and validation of further 
approximate treatments for example in the spirit 
of (TD)OEP approaches. A first application to 
a one-dimensional molecule proves the feasibility 
and stability of the scheme. A comparison with 
exact exchange confirms the need and satisfying 
performance of full TDSIC in ionization dynam- 
ics. 
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We work in the Kohn-Sham scheme of DFT, buht 
on a set of single particles wavefunctions {ipa, a: = 
1,...,A^}. In the SIC case, the starting energy 
functional reads 
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where the electronic kinetic energy E'kin is com- 
plemented by the external ionic contribution -Eion 
and the LDA approximation -Elda [p] to the exact 
electron-electron interaction energy. The densities 
are further defined as pa = IV'qP and p = pa- 
Note that all summations run over occupied states 
only. 

The TDSIC equations are obtained from the prin- 
ciple of stationary action using the SIC energy 
functional ([1]) 
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within explicitely including the orthonormality 
constraint with Lagrange multipliers Vari- 
ation of i^sic with respect to V'q leads to single 
particle equations for the in which the one- 
body Hamiltonian is obtained as 
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hatpa , ha = ft-LDA ~ Ua , (3a) 
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The emerging one-body Hamiltonian depends 
on the state ipa on which it acts through the SIC 
term Uq.. The /i^'s can be simply recast in a SIC 
Hamiltonian /isic by employing projectors as 



'iSIC = ^LDA " t/a|'0a)(V'o 
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That form embodies the state-dependence in the 
projector and displays clearly the non-hermitian 
nature of /isiC: which is also not invariant under a 
unitary transform amongst the iV'a)- 



Variation of the action with respect to t/"* thus 
yields the TDSIC equations as 

{hsic - ifidt) IV'q) = J2p IV'/3)^/3a : (4a) 

A/3q = (-0/3 1 - i^^tlV'a) ■ (4b) 

together with the symmetry condition 

0=i4'f3\Up-Ua\4'c.). (4c) 

which has to be fulfilled at each instant. It should 
be noted that once one has achieved the symmetry 
condition (|lc)) . the SIC Hamiltonian ([5c)l acquires 
an interesting property. Although it remains non- 
hermitian as a whole, it becomes hermitian within 
the space of occupied states 



(^l/lSIclV'a) = i:>P alhsicl^pY 
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The above TDSIC equations are quite involved as 
the time propagation, Eqs. (|4all4bp . is constrained 
by the symmetry condition (|^c)) . unlike propaga- 
tion with a strictly hermitian Hamiltonian (LDA 
or Hartrce-Fock) for which the symmetry condi- 
tion is automatically fulfilled. A simple minded 
step 

|V'a(<)) =CXp(-^ / dt'/lsic(i')Wa(<0)), 



is thus not directly applicable because it does not 
ensure the preservation of the symmetry condition. 
One has to employ an interlaced step which allows 
to fulfill simultaneously Eqs. (jH). 

In order to overcome this difficulty, we note that 
there is always the freedom of unitary trans- 
formations amongst the set of occupied orbitals 
{tpa,ct = l,...,iV} without changing the state of 
a system. The ipa are the unique ingredients of 
the SIC mean field. But propagation of the whole 
state may be formulated in terms of another set 
{(fi} chosen to have convenient propagation prop- 
erties and connected to the ipa by a unitary trans- 
formation within occupied states : 
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We call the {(pi} the "propagating" set and the 
{V'q} the "symmetrizing" set, the idea being to 
perform a joined propagation of {^i} and {ipa} 
such that each of the two (connected) sets of or- 
bitals contributes either the propagation or the 
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symmetry condition. It is nevertheless crucial to 
note that the V'a remain the key constituents com- 
posing the SIC potentials ([5cl) and entering the 
symmetry condition (|4c|) through the Ua- The 
propagation set ipi gives the freedom to choose 
the propagation within occupied states at conve- 
nience. We choose it such that 



(6b) 



which allows to fulfill Eqs. (|4all4bp and which is 
possible as soon as the V'q fulfill the symmetry 
condition (Pc|) . This allows then to propagate the 
{(Pi} as : 



/isic - ifidt ]\ipi) = 
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and to care for the symmetry condition at the side 
of the tpa, the latter fixing the coefficients of the 
unitary transformation. This reads formally 
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It is to be noted that the propagator in Eq. ([6c)) 
is not strictly unitary because /isic is not hermi- 
tian. But the hermiticity within occupied space, 
Eq. © , guarantees that the propagation ([5c|) pre- 
serves orthonormality within occupied space, i.e. 
{ip,{t)\tpj{t)) = S,j. 

The above described propagation scheme for TD- 
SIC exhibits some interesting properties. As al- 
ready noted, it preserves orthonormality, which 
is a crucial requirement for any time evolution 
scheme. One can also show that energy and the 
zero-force theorem are conserved (as long as there 
are no time-dependent external fields). Finally, we 
remark that the choice ([6cl) for the propagating set 
is not the only possibility. There is some freedom 
for other choices as, e.g., optimizing single-particle 

94 

energies". 

As a final point, we apply TDSIC to a typical 
example of laser excitation and subsequent ion- 
ization of a molecule. We use a simple one- 
dimensional model for a molecule in the spirit 
o:^. The model case is a dimer with two elec- 
trons having the same spin. As an interaction, 
we use the following smoothed Coulomb potential, 

Wij — I \J {x\ — + a?j , where the parameters 
flij for electron-electron, electron-ion and ion-ion 
interactions are tuned to reproduce typical molec- 
ular energies. Taking that interaction, we develop 
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FIG. 1: Time evolution of the number of emitted elec- 
trons for a dimer molecule with two electrons. Com- 
pared are results from full TDHF, TDLDA and TD- 
SIC. 



with LDA an energy functional for the exchange 
term. Working at the level of exchange only allows 
to have fully fiedged time-dependent HE (TDHF) 
calculations as benchmark to which DET calcula- 
tions can be compared. 

A very short laser pulse is simulated as an instan- 
taneous boosts^. This has the advantage that en- 
ergy conservation can be used as test for the calcu- 
lations. We have checked conservation of energy, 
orthonormality and symmetry condition which are 
all found fully satisfied. As a further observable, 
we consider the degree of ionization which, as 
stated above, is a sensitive quantity to probe the 
effect of SIC. The results are shown in Fig. [1] com- 
paring the TDHF benchmark with TDLDA and 
TDSIC. It is obvious that TDSIC comes very close 
to the benchmark. 

Finally we want to remark that the double-set 
strategy can also be applied to the static SIC 
problem. The static SIC equations are obtained 
by minimizing Ssic with a constraint on or- 
thonormality of the single-particle wavefunctions 

■00, {Esic - Y.i3,^{i'0\'^i)^ip) = 0' following 
standard variational derivation a^^'^'^ . One then 
simply obtains /isiclV'a) = Z]/3 \i'i3)^ap, again to- 
gether with the symmetry condition Ea. (|4c)) . Note 
that the matrix of Lagrange multipliers X^a = 

(^/3|/isic|'0q) is usually non-diagonal. The states 
ipa which emerge as solutions of stationary equa- 
tions are optimized to produce the correct SIC 
potentials. They do not give any clue on single- 
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particle energies. One can now introduce a second 
set of {ipi} connected by the stationary analogue of 
the unitary transformation (j6ap and require, e.g., 
that these ipi diagonalize ft-sic (expressed in terms 
of the ipa, Eg. iPc]) ). or equivalently the matrix of 
Lagrange multipliers Xafj. The eigenvalues thus 
obtained can be interpreted as single-particle en- 
ergies and they are found to agree fairly well with 
HF values in the case of exchange only calcula- 
tions that we consider here. Note finally that these 
static SIC equations also emerge naturally as the 
stationary limit of TDSIC. 

We have proposed in this letter a consistent vari- 
ational formulation of time-dependent SIC (TD- 
SIC) together with a manageable and transparent 
scheme for the solution of the TDSIC equations. 
This scheme provides conservation of energy, zero- 
force theorem, and orthonormality of the occupied 
states. The stationary limit of static SIC is also 
properly recovered both by the theory itself and by 
the propagation scheme. We applied the scheme 
to laser-induced ionization of a one-dimensional 
dimer molecule as a test case. The calculations 
have proven to run stable and to fulfill all theoret- 
ical constraints. As a critical observable, we have 



investigated the time evolution of ionization and 
found nice agreement of the exchange only TDSIC 
results with an exact TDHF calculation. 

Although even full 3D calculations have proven 
to be feasible, it is to be admitted that fully 
fledged TDSIC is rather involved and thus com- 
putationally complex. We consider it as start- 
ing point for further development towards sim- 
plified schemes. A promising option is provided 
by a time-dependent form of optimized effec- 
tive potential o^^'^° . Work in that direction is in 
progress. Nevertheless the full TDSIC serves as a 
crucial benchmark for such developments. 
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